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Abstract The effects of viscosity on the circumplanetary disks residing in the vicinity of 
protoplanets are investigated through two-dimensional hydrodynamical simulations with 
■ the shearing sheet model. We find that viscosity can affect properties of the circumplane- 

tary disk considerably when the mass of the protoplanet is Alp < SSM®, where is the 
Earth mass. However, effects of viscosity on the circumplanetary disk are negligibly small 
when the mass of the protoplanet Mp > 33 Af^. We find that when Mp < SSA/q, vis- 
cosity can disrupt the spiral structure of the gas around the planet considerably and make 
the gas smoothly distributed, which makes the torques exerted on the protoplanet weaker 
Thus, viscosity can make the migration speed of a protoplanet lower After including vis- 
cosity, size of the circumplanetary disk can be decreased by a factor of > 20%. Viscosity 
helps to transport gas into the circumplanetary disk from the differentially rotating cir- 
cumstellar disk. The mass of the circumplanetary disk can be increased by a factor of 
50% after viscosity is taken into account when Alp < 33A/0. Effects of viscosity on the 
formation of planets and satellites are briefly discussed. 

Key words: accretion, accretion disks — hydrodynamics — planets and satellites; for- 
mation — Solar system: formation 

1 INTRODUCTION 

Up to date, more than 500 exoplanets have been detected. Most of the exoplanets are gas giant planets, 
as massive planets are preferentially observed by current detection methods. Thus, it is important to 
understand the formation process of gas giant planets. According to the core accretion model, a solid 
core with several A/®, forms first through coagulation of planetesimals in the circumstellar disk. The 
protoplanet captures a hydrostatic envelope when its mass is less than Alp < IOA/q (Mizuno 1980; 
Stevenson 1982; Bodenheimer & Pollack 1986; Pollack et al. 1996; Ikoma, Nakazawa & Emori 2000; 
Ikoma, Emori & Nakazawa 2001 ; Hubickyj, Bodenheimer & Lissauer 2005). Ikoma et al. (2000) showed 
that run-away gas accretion is triggered when the solid core mass exceeds ~ 5 — 20A{[q, the protoplanet 
quickly increases its mass by gas accretion. A gas giant planet acquires almost all of its mass in the 
run-away gas accretion phase. 

Since gas accreting from a differentially rotating circumstellar disk has nonzero angular momentum, 
a circumplanetary disk can form around the protoplanet (Tanigawa & Watanabe 2002). The circumplan- 
etary disk can influence many aspects of the protoplanet. For example, previously, when calculating 
the torque on the protoplanet, the contribution of gas inside the whole Roche lobe (or Hill radius) is 
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neglected. This may be not appropriate. The size of the circumplanetary disk may be smaller than the 
Roche lobe (as shown in Section 3.2); therefore, when calculating the torque on the protoplanet the gas 
inside the Roche lobe but beyond the outer boundary of the circumplanetary disk should be taken into 
account. Also, properties of the circumplanetary disk can determine the evolution of protoplanet and its 
resulting mass. Finally, satellites form in the vicinity of the protoplanet; studying the properties of the 
circumplanetary disk helps to investigate the formation process of the satellites. 

Gas accretion process onto a protoplanet have been investigated in global simulations by many 
authors (e.g. Bryden et al. 1999; Kley 1999; Lubow, Seibert & Artymowicz 1999; Kley, D'Angelo & 
Henning 2001; D'Angelo, Henning & Kley 2002; Bate et al. 2003; D'Angelo, Kley & Henning 2003). 
However, since the main purpose of these studies was gap formation and planet migration on a large 
scale, the region in the vicinity of the protoplanet has not been investigated with sufficient resolution. 
Thus, in their simulations, properties of the circumplanetary disk were not thoroughly explored. The 
fine structure of the circumplanetary disk has been investigated with shearing sheet models without 
viscosity (Tanigawa & Watanabe 2002; Machida et al. 2008; Machida 2009; Machida et al. 2010). A 
question is the mechanism of angular momentum transport in the circumplanetary disk in their models. 
Gravitational interaction between the protoplanet and gas can produce spiral shocks inside the Roche 
lobe (or Hill radius) of a protoplanet. Gas flows into the Hill sphere of a protoplanet through the inner 
and outer Lagrange points. The gas that flows into the Roche lobe from the inner (outer) Lagrange point 
will undergo a strong shock on the opposite outer (inner) Lagrange side of the Roche lobe, angular 
momentum is lost through the collision between gas and shocks. The gas spirals inward toward the 
protoplanet as a result of successive shocks. 

Obviously, there should be viscosity in the circumstellar disk, which drives the gas flow in the 
disk toward the central star The most promising origin of viscosity in the circumstellar disk should 
be magnetic turbulence generated by the magnetorotational instability (MRI) (Balbus & Hawley 1991; 
1998). Previous work found that despite the ionization rate of a circumstellar disk is low, magnetic 
field can remain dynamically important, and MRI perturbations can grow under a wide range of fluid 
conditions and magnetic field strengths (Salmeron & Wardle 2005). Viscosity may play important roles 
on the properties of the circumplanetary disk. The angular momentum profiles of a circumplanetary disk 
may be affected significantly by viscosity. Also, viscosity helps to transport gas into the circumplanetary 
disk from the differentially rotating circumstellar disk, the mass of the circumplanetary disk may be 
affected by viscosity significantly. Therefore, it is of great importance to study the effects of viscosity 
on properties of circumplanetary disk. 

In this paper, we study the effects of viscosity on circumplanetary disks with the shearing sheet 
models. We use an anomalous stress tensor to mimic the shear stress originated from magnetohydrody- 
namic (MHD) turbulence. In Section 2, we describe our models. Results are described in Section 3. We 
discuss and summarize our results in Section 4. 



2 MODEL 
2.1 Equations 

We confine our models to two-dimensions. We assume that the temperature is constant and that the self- 
gravity of the disk is negligible. The orbit of the protoplanet is assumed to be circular in the equatorial 
plane of the circumstellar disk. The protoplanet is not allowed to migrate. 

We consider a local region around a protoplanet, using the shearing sheet model (e.g. Goldreich & 
Lynden-Bell 1965). We take local Cartesian coordinates rotating with the protoplanet with the origin at 
the protoplanet and the x~ and y— axes in the radial and azimuthal direction of the disk, respectively. 
We solve the equations of hydrodynamics without self-gravity: 

as , , 

— + V-(Sv)=0 (1) 
9v 1 1 

— + (v • V)v = --VP - V$ - 2r!pe, x v + -V ■ T (2) 
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where E is the surface density, v is the velocity, P is the vertically integrated gas pressure, $ is the 
gravitational potential, flp is the Keplerian angular velocity of the protoplanet, Cz is the unit vector 
along the rotation axis of the protoplanet, T is the vertically integrated anomalous stress tensor We 
adopt an isothermal equation of state, 

P = (3) 
where Cs is the sound speed. The angular velocity of the protoplanet is given by 

/GM, ^ 



^P={-j^] (4) 



where G, Mc and Op are gravitational constant, the mass of the central star and the orbital radius of the 
protoplanet, respectively. The gravitational potential is given by 

$ = --Wy ^ (5) 

2 r 

where Mp and r are the mass of the protoplanet and the distance from the center of the protoplanet, re- 
spectively. The first term is composed of the gravitational potential of the central star and the centrifugal 
potential. The second term is the gravitational potential of the protoplanet. The Hill radius inside which 
the protoplanet gravity dominates is defined as 



Using the Hill radius, we can rewrite equation (4) as 

We use the stress tensor T to mimic the shear stress, which is in reality magnetic stress associated 
with magnetohydrodynamics (MHD) turbulence driven by the magnetorotational instability (MRl). We 
assume that the only non-zero component of the stress tensor T is 

This is because MRl is driven only by the shear associated with orbital dynamics. In equation (7), 
ji ~ Sac^/rij.. a, Cs and $7^ are viscosity coefficient, sound speed and Keplerian angular velocity about 
the protoplanet, respectively. Previous paper (e.g. Papaloizou, Nelson & Snellgrove 2004) studying the 
interaction between a protoplanet and magnetized circumstellar disk found that a ~ 3 x lO^'^. The 
maximum a used in this paper is 3 x 10^"^. 

In this paper, we normalize length by the scale height of the circumstellar disk h = Cs/^p, time 
by the inverse of the Kplerian angular velocity of the protoplanet and the surface density by the 
unperturbed surface density of the standard solar nebular model (Hayashi 1981; Hayashi et al. 1985). 



2.2 Numerical method 

The numerical simulations are performed using the ZEUS-2D code (Stone & Normal 1992a; 1992b). 
Our initial settings are similar to those of Machida et al. (2008). The gas flow has a constant shear in the 
x— direction as 




Initially, the gas has uniform surface density of the unperturbed disk. In this paper, our standard com- 
putational domain is that \x\ < 6(= Xmax) and \y\ < 12(= Umax)- We adopt logarithmic spacing grids 
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with the finest resolution {Ax = 0.003) around the protoplanet. The total number of grids for our stan- 
dard computational domain is 500 x 1000, the resolution is high enough to study the circumplanetary 
disk in the vicinity of the protoplanet. We inject gas with the linearized Keplerian shear on ?/ = Vmax, 
{0 < X < Xmax) and y = —ymax, {—Xmax < X < 0). Foi the rest of the boundaries, we adopt outflow 
boundary conditions. 

In order to avoid singularity in the proximity of the protoplanet, the gravitational potential of the 
protoplanet is smoothed in its neighborhood with 

where rsm is the smoothing length of the protoplanet's potential. In this paper, we choose rgm = 0.05. 
Tanigawa & Watanabe (2002) showed that rgm = 0.05 is safe to study the properties of the circumplan- 
etary disk in the vicinity of the protoplanet. 

In the run-away gas accretion phase, accretion by planet is non-negligible. There should be accre- 
tion, neglecting accretion is unphysical. Therefore, we should mimic the accretion process by planet. 
Because the growth timescale of a planet is much longer than the typical time of our simulations, we 
should have a constant accretion rate which is independent from the parameters used to mimic the ac- 
cretion process. As done by Tanigawa & Watanabe (2002), the gas inside rgmk is removed by a constant 
rate [1]"+^ = S"(l — At)], At(^ 1) is a time step of the calculations, and superscript n is the number 
of numerical time step. Tanigawa & Watanabe (2002) have tested the effects of r^m and rsink on the 
accretion rate of the planet. They found that the results do not depend on the values of Vsm and rsink as 
long as rs,n = rsink < 0.07. In fliis paper, we set rsm. = rsink = 0.05. 

2.3 scaling 

In the standard solar nebular model (Hayashi 1981; Hayashi et al. 1 985), the temperature T, sound speed 
Cs, and gas density po are given by 

( lV'^ ( a„ \-i/2 
r = 280 — (— (11) 



Lq) \\AV I 
where L and Lq are the protostellar and solar luminosities; 

where [i = 2.34 is the mean molecular weight of the gas composed mainly of Hi and iJe; and 

PO = 1.4 X 10-9 (^) g/cm^, (13) 

respectively. When the values of = IAIq and L = ILq are adopted, using equations (4), (11) and 
(12), we can describe the scale height h as 

/a \ 5/4 

/i = 5.0 X 10" ( — ^ cm. (14) 
\1AUJ 

The mass of the planet in unit of Jupiter mass can be described as 

-o.uf^y'^U^Y'U^X (15) 



Mjup \IMqJ \1AU/ V h , 

In this paper, we assume that Mc — IMq, L — ILq. The planet is located at Op = 5.2AU. 
Therefore, the temperature of the gas is T = 123K. Because our shearing box just represent a local 
region around the planet, we further assume that the temperature is uniform in the whole computational 
domain. 

In the paper below, the Hill radius Rh is in unit of h. In table 2, Rh = 0.4, 0.5 and 0.63 correspond 
to planet mass of 0.026, 0.05 and 0.1 Jupiter mass, respectively. 
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Fig. 1 Specific angular momentum (upper-panel) and 
surface density (lower-panel) of the circumplanetary 
disks for a IG.SA/g protoplanet as a function of dis- 
tance r from the protoplanet with different sheet size. 
It can be seen that the effects of sheet size on the struc- 
ture of the circumplanetary disk are negligible as long 
as the sheet size is much larger than the Hill radius. 



Table 1 Parameters for the test models 



Models 


Sheet size 


Grids 


Rh 


a 


MOOSVO 


12 X 24h 


500 X 1000 


0.5 





BM005V0 


24 X iSh 


656 X 1312 


0.5 





SM005V0 


6 X 12h 


358 X 716 


0.5 






2.4 Tests of the effects of sheet size on the results 



In this section, we study the effects of sheet size on the structure of the circumplanetary disk inside 
the Hill radius. We find that the flow structure inside the Hill sphere (or circumplanetary disk) does not 
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Fig. 2 Contours of logarithm surface density of the 
circumplanetary disks surrounding protoplanets with 
mass of 8.5, 16.5 and SSM^ (top to bottom rows). 
The top row panels show the circumplanetary disks 
surrounding a S.SM® planet. In the top row, the left- 
, middle- and right-hand panels correspond to mod- 
els M0026V0 {a = 0), M0026V1 {a = lO^^) and 
M0026V2 (a = 3 X lO'^), respectively. The mid- 
dle row panels show the circumplanetary disks sur- 
rounding a 16.5Af0 planet. In the middle row, the left- 
, middle- and right-hand panels correspond to models 
M005V0 (a = 0), M005V1 {a = lO^^) and M005V2 
(a = 3 X lO^'^), respectively. The bottom row panels 
show the circumplanetary disks surrounding a 337\/0 
planet. In the bottom row, the left-, middle- and right- 
hand panels correspond to models MOIVO (a = 0), 
MOIVI (a = 10-3) and M01V2 (a = 3 x IQ-^), 
respectively. 



depend on the computational box size as long as the size is much larger than the Hill radius. Here we 
just take Rh = 0.5, a = as a example. Table 1 lists the parameters of the Rh = 0.5 tests. 

Fig.l shows the specific angular momentum (upper-panel) and surface density (lower-panel) distri- 
bution of the circumplanetary disks for a 16.5A/0 protoplanet with different sheet size. From this figure, 
we can see that the effects of sheet size on the structure of the circumplanetary disk are negligible as 
long as the sheet size is much larger than the Hill radius. The results shown in Section 3 is calculated 
with our standard sheet size 12 x 24/i. 
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Table 2 Parameters for all of our models 



Models 


Rh 


Mp 




a 




M0026V0 


0.4 


8.5Me 









M0026V1 


0.4 


8.5Me 




10-3 




M0026V2 


0.4 


8.5Me 


3 


X IQ- 


■3 


M005V0 


0.5 


16.5Me 









MOOS VI 


0.5 


16.5Me 




10-3 




M005V2 


0.5 


16.5Me 


3 


X 10" 


-3 


MOIVO 


0.63 


33M© 









MOIVI 


0.63 


33M0 




10-3 




M01V2 


0.63 


33Me 


3 


X 10~ 


3 



3 RESULTS 

Table 2 lists all of the models in this paper. We find that all of the models have settled into their steady 
state before 20 orbits, so that the results shown below at i = 75 orbits are fully relaxed. 

In order to investigate the effects of viscosity on the circumplanetary disk, for each protoplanet with 
different mass, we carry out three simulations with different viscosity coefficient a. We find that when 
the protoplanet mass Mp > SSAf^, the effects of viscosity are negligible. 

3.1 Circumplanetary disk structure and migration of protoplanets 

Fig. 2 shows the circumplanetary disks structure around protoplanets with mass of 8.5, 16.5 and SSA/^ 
(top to bottom rows). The top row panels show the circumplanetary disks surrounding a 8.5A/0 proto- 
planet. In the top row, the left-, middle- and right-hand panels correspond to models M0026 VO {a = 0), 
M0026V1 (a = 10"^) and M0026V2 (a = 3 x 10"^), respectively. The middle row panels show 
the circumplanetary disks surrounding a 16. SM® protoplanet. In the middle row, the left-, middle- 
and right-hand panels correspond to models M005V0 (a = 0), MOOS VI (a = IQ-^) and M005V2 
(a = 3 X 10^"^), respectively. The bottom row panels show the circumplanetary disks surrounding a 
33A/0 protoplanet. In the bottom row, the left-, middle- and right-hand panels correspond to models 
MOIVO {a ^ 0), MOl VI (a = IQ-^) and M01V2 (a 3 x IQ-^), respectively. It can be seen clearly 
that as the mass of the protoplanet increases, the spiral structure of the circumplanetary disk becomes 
more prominent, higher mass protoplanet can excite higher amplitude spiral shock. 

For the protoplanet with mass of 8.5A/©, with the increase of strength of viscosity, the spiral struc- 
ture is getting weaker The spiral structure completely disappears when q = 3 x 10"^ and the circum- 
planetary disk is nearly axisymmetric about the protoplanet. For a 8.5A/0 protoplanet, the maximum 
surface density of the spiral waves is only 2 times bigger than the minimum surface density; the am- 
plitude of the spiral wave is weak, viscosity can easily disrupt the spiral structure and make the gas in 
the circumplanetary disk smoothly distributed. For the protoplanet with mass of 16. SA/^, viscosity also 
makes the spiral structure of the circumplanetary disk weaker But the effect is small compared to the 
smaller protoplanet mass case (Alp = 8.5 A/®). This is because with the increase of the mass of the pro- 
toplanet, the amplitude of the spiral waves gets higher, viscosity can hardly affect the spiral structure. 
It can be seen that when the protoplanet mass reaches 33A/0, viscosity almost plays no roles on the 
circumplanetary disk structure. 

The protoplanet excites spiral density waves at the Lindblad resonance and the torques exerted on 
the protoplanet as the reaction of exciting waves make the orbit of the protoplanet around central star 
change. Previous works (Lubow et al. 1999; Tanigawa & Watanabe 2002) found that the spiral structure 
around the protoplanet can affect the torques significantly. In the limiting case, if the gas distribution 
around a protoplanet is perfectly axisymmetric about the protoplanet, the torques exerted on the pro- 
toplanet will be zero. From Fig. 2, we find that viscosity can disrupt the spiral structure considerably 
when the mass of protoplanet Mp < 33A/0. Thus, we can expect that when the protoplanet is small 
{Mp < 33A/0), viscosity can affect the torques exerted on the protoplanet considerably. 
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It is always believed that the gas inside the Hill sphere of a planet migrates with the planet. 
Therefore, in most previous works, when calculating the torque exerted on a planet, the contribution 
from the gas inside the Hill sphere is completely excluded. However, Crida et al. (2009) showed 
that the gas bound to the planet (circumplanetary disk) only exists inside Q.5Rh- In the region 
0.5Rh < r < Rff, the gas is not bound to the planet. Therefore, when calculating the torque, the 
gas inside the Hill sphere should not be excluded completely. 

Now, we quantitatively study the effects of viscosity on the torques exerted on a protoplanet. Since 
we adopt the local approximation, the net torque in our simulation becomes exactly zero. But an actual 
net torque is not zero, owing to slight asymmetry of density distribution, temperature and the curvature 
of the protoplanet orbit. The net torque exerted on the protoplanet is roughly proportional to the one- 
side torque exerted by the gas exterior to the orbit of the protoplanet (Tanigawa & Watanbe 200213; 
therefore, the one-side torque would be a clue to solve the migration problem. Because we adopt the local 
approximation, we just discuss the effects of viscosity on the one-side 'force' exerted on the protoplanet. 
The 'force' corresponds to the torque divided by the semimajor axis of the protoplanet. The one-side 
force exerted on a protoplanet is defined by 

Fy^out{rmask) = / / T:—0{r - r,nask)dydx (16) 

where r„iask is the artificial inner limit of the force integration, and 6 is the step function. If r > rmask, 
9=1, otherwise 9 = 0. The 'force density' is defined as dF{x) = J^™"'" Y,^dy. 









^ -0 

o 



-0 
-0 

0.1 1.0 

r 

mask 

Fig. 3 Gravitational force exerted on the 8.5M© pro- 
toplanet from disk gas in a; > as a function of 
"Tmask, within which gas is excluded from force inte- 
gration. In this figure, the solid, dotted and dashed lines 
correspond to models M0026V0 (a = 0), M0026V1 
(a = 10^3) jjjj^j M0026V2 (a 3 x IQ-^), respec- 
tively. 

We show Fy^out as a function of Vmask for a protoplanet with mass of S.SM® in Fig. 3. In Section 
3.2, we show that the outer boundary of S.SAf® and 16.57\/0 protoplanets located inside r = 0.22, 
therefore we can focus on the torque exerted by the gas beyond r = 0.22. The gas in a; > exerts a 




' Net torque is the difference of the two opposite-signed one-side torques, and the ratio of the two one-side torques does not 
change very much when one change parameters such as planet mass or temperature in linear regime. Thus, at least in this regime, 
we can say that the net torque is proportional to one-side torque. 
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X 

Fig. 4 Upper-panel shows the gravitational force ex- 
erted on the 16.5Af0 protoplanet from disk gas in 
x > as a function of r„iask, within which gas is ex- 
cluded from force integration. In this panel, the solid, 
dotted and dashed lines correspond to models MOOS VO 
(a = 0), M005V1 {a = lO^^) and M005V2 (a = 
3 X 10""^), respectively. Lower-panel shows the torque 
density as a function of distance from the 16.5A/0 pro- 
toplanet. In this panel, the solid, dotted and dashed 
Unes correspond to models M005V0(q; = 0), MOOS VI 
(a = 10^3) and MOOS V2 (a = 3 x lO"^) respectively. 

negative torque on the protoplanet, which pushes the protoplanet migrating towards the central star. As 
expected, when the spiral structure around the planet gets weaker, the force exerted on the protoplanet 
is decreased significantly. Fig. 4 (upper-panel) plots Fy^out as a function of Vmask for a protoplanet with 
mass of 16.5Af0. When a = 3 x 10^"^, the spiral structure is partially disrupted by viscosity, which also 
leads to the reduction of the torque exerted on the planet. In order to show clearly the 'torque' is affected 
by viscosity, the lower panel of Fig.4 plots the 'torque density' as a function of distance from the planet. 
The case for a S.SM^ protoplanet is similar. We find when Mp > 33Af0, the effect of viscosity on 
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the torque exerted on protoplanets is negligible. We conclude that the viscosity can affect the torques 
exerted on a protoplanets considerably by disrupting the spiral structure of gas disk when the mass of 
the protoplanet is small [Mp < SSA/^). 



1 .00 F 




r 



Fig. 5 Specific angular momentum of the circumplan- 
etary disks surrounding a 8.5A/© protoplanet. In this 
figure, the solid, dotted and dashed lines correspond to 
models M0026V0 {a = 0), M0026V1 (a = IQ-^) 
and M0026V2 (a = 3 x 10"^), respectively. The long- 
dashed line corresponds to the Keplerian angular mo- 
mentum with respect to the protoplanet. 



3.2 Circumplanetary disk size 

Now, we investigate the effects of viscosity on the size of the circumplanetary disks. The radial edge of 
a circumplanetary disk is taken as the point of turnover in the specific angular momentum of the disk. 
Quillen & Trilling (1998) made a simple analytic prediction of the approximate circumplanetary disk 
radii around accreting protoplanets. They assume that the gas flows into the Hill sphere of a protoplanet 
via the inner and outer Lagrange points. They also assume that when the gas arrives at the Lagrange 
points, the velocity, relative to the Lagrange points, is negligibly small. The Lagrange points have the 
same angular velocity around the central star as the protoplanet. The Lagrange points are approximately 
located at 

r = ap±RH (17) 

Thus, when the gas at the Lagrange points is captured by the protoplanet, its specific angular mo- 
mentum relative to the protoplanet is 

J - Rl^p (18) 

Assuming conservation of angular momentum when the accreting gas flows towards the protoplanet, the 
centrifugal radius, r^, of the circumplanetary disk is 

- ^ (19) 



which yields 



Rnli- (20) 
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Fig. 6 Specific angular momentum of the circumplan- 
etary disks surrounding a 16.5M^ protoplanet. In this 
figure, the solid, dotted and dashed lines correspond 
to models M005V0 (a = 0), MOOS VI (a = IQ-^) 
and M005V2 (a = 3 x 10"^), respectively. The long- 
dashed line corresponds to the Keplerian angular mo- 
mentum with respect to the protoplanet. 



Thus, the typical size of a circumplanetary disk should be approximate Rh /3. 

Fig. 5 shows the specific angular momentum of the circumplanetary disks around a S.SMg pro- 
toplanet. In this figure, the solid, dotted and dashed lines correspond to models M0026V0 {a = 0), 
M0026V1 (a = 10"^) and M0026V2 (a = 3 x 10^^), respectively. We can see the outer boundary of a 
non-viscous circumplanetary disk (model M0026V0) around a 8.5Mq protoplanet locates at r = 0.15, 
the disk size does not differ greatly from the Rh /3 prediction. Beyond r = 0.15, the gas rotates around 
the central star, the angular momentum becomes smaller and smaller with increasing distance from the 
protoplanet. The angular momentum relative to the protoplanet even becomes negative when the dis- 
tance from the protoplanet is sufficiently big which is not shown in Fig. 5. When viscosity is included, 
the size of the circumplanetary disk decreases. In model M0026V2, the outer boundary of the circum- 
planetary disk locates at r ~ 0.11, the size of the disk is reduced by a factor of 27% compared to the 
non-viscous disk. As can be seen from Fig. 5, the specific angular momentum of the circumplanetary 
disk becomes lower after including viscosity, which is due to the outward angular momentum transfer 
by viscosity. The variation of angular momentum at given radius is determined by the divergence of the 
viscous stress; at the outer boundary of the circumplanetary disk, the divergence of viscosity is stronger 
than that at other radii due to the stronger shear of the gas (it can be clearly seen from Fig. 5), which 
makes the decreasing amplitude of the angular momentum at the outer boundary much larger than that 
of the gas at other radii. Therefore, the outer boundary of the circumplanetary disk is moved inward 
when viscosity is taken into account. 

Fig. 6 shows the specific angular momentum of the circumplanetary disks around a 16.5Af0 proto- 
planet. In this figure, the solid, dotted and dashed lines correspond to models M005V0 (a = 0), MOOS VI 
(a ~ 10^^) and MOOS V2 (a = 3 x 10^"^), respectively. The outer boundary of the circumplanetary disk 
is also moved inward when viscosity is included. The outer boundary of a non-viscous circumplanetary 
disk (model MOOSVO) around a 16.5M0 protoplanet locates at r = 0.22. In model M00SV2, the outer 
boundary of the circumplanetary disk locates at r = 0.17, the size of the disk is reduced by a factor of 
23% compared to the non-viscous disk. We find that when the protoplanet mass Alp > 33M0, the effect 
of viscosity on the size of a circumplanetary disk is negligibly small. 
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Fig. 7 The surface density averaged over y against x 
for the gas around protoplanets. The upper and lower 
panels correspond to 8.5 and 16.5M0 protoplanets, re- 
spectively. In the upper panel, the solid, dotted and 
dashed lines correspond to models M0026V0 (a = 0), 
M0026V1 (a = IQ-^) and M0026V2 (a = 3 x IQ-^), 
respectively. In the lower panel, the solid, dotted and 
dashed lines correspond to models M005V0 {a = 0), 
MOOS VI {a = 10-3) and M005V2 (a = 3 x IQ-^), 
respectively. 

Crida et al. (2009) have found that the size of the circumplanetary disk ~ 0.5Rh, which is slightly 
larger than that obtained in this paper. This is because Crida use a energy equation to solve internal 
energy of the gas. With the energy equation, the collapse of the gas onto the circumplanetary disk is 
limited by the heating due to adiabatic compression, which gives a wider circumplanetary disk than in 
the locally isothermal case. 

From Fig. 5 and Fig. 6, we see that the rotational velocity of the disk around the planets is signifi- 
cantly sub-Keplerian. The reason is as follows. In the radial direction, gravitational force is balanced by 
centrifugal force and pressure gradient force. If the temperature of the gas is fixed, with the increasing 
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Fig. 8 Surface density of the circumplanetary disks 
as a function of distance r from the protoplanet. 
The upper and lower panels correspond to 8.5 and 
I6.5M0 protoplanets, respectively. In the upper panel, 
the solid, dotted and dashed lines correspond to mod- 
els M0026V0 (a = 0), M0026V1 {a IQ-^) and 
M0026V2 (a = 3 X 10"^), respectively. In the lower 
panel, the solid, dotted and dashed lines correspond to 
models M005V0 (a = 0), M005V1 (a = lO^^) and 
M005V2 (q' = 3 X 10^3), respectively. 



of planet mass, the centrifugal force will increase (rotational velocity increase). In our paper, the mass 
of the planet is small (around 10 earth mass), so pressure gradient force is important, the gas rotates 
sub-keplerian. Tanigawa, Ohtsuki & Machida (2012) found that the gas around the planet rotates almost 
Keplerian. The reason is as follows. In Tanigawa, Ohtsuki & Machida, ApJ (2012), the temperature is 
identical to that in our paper, but their planet mass is 1 Jupiter mass which is much bigger than the planet 
mass in our paper Therefore, in Tanigawa, Ohtsuki & Machida (2012), the gas rotates almost keplerian. 
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3.3 Circumplanetary disk mass 

In the standard picture, a protoplanet residing in a circumstellar disk can exert torques on the circumstel- 
lar disk through the excitation of spiral density waves (e.g. Goldreich & Tremaine 1979). The angular 
momentum carried by the density waves will be deposited in the circumstellar disk where the waves are 
damped. The circumstellar disk gas exterior to the protoplanet orbit gets positive angular momentum 
and thus moves outward. The circumstellar disk gas interior to the protoplanet orbit gets negative angu- 
lar momentum and thus moves inward. For a low-mass protoplanet, a partial density gap forms along 
the orbit of the protoplanet. However, the viscosity inside the circumstellar disk tries to transport gas 
into the low density partial gap region, which makes the density gap less prominent. Viscosity can affect 
the depth of the density gap around the orbit of the protoplanet and thus affects the circumplanetary disk 
mass. 

In Fig. 7, we plot the surface density averaged over y against x for the gas around protoplanets. The 
upper and lower panels correspond to 8.5 and IG.SAf® protoplanets, respectively. In the upper panel, 
the solid, dotted and dashed lines correspond to models M0026V0 (a = 0), M0026V1 (a = 10"^) and 
M0026V2 (a = 3 X 10^^), respectively. In the lower panel, the solid, dotted and dashed lines correspond 
to models M005V0 (a = 0), MOOS VI (a = IQ-^) and M005V2 (a = 3 x IQ-^), respectively. It is 
obviously seen that the protoplanets try to open a density gap along its orbit. As expected, the protoplanet 
in the inviscid gas produces a deeper and wider gap than others. 

In Fig. 8, we plot the surface density profiles of the circumplanetary disks around protoplanets. 
The upper and lower panels correspond to 8.5 and 16.5Af® protoplanets, respectively. In the upper 
panel, the solid, dotted and dashed lines correspond to models M0026V0 {a = 0), M0026V1 {a = 
10"'^) and M0026V2 (a = 3 x 10^"^), respectively. In the lower panel, the solid, dotted and dashed 
lines correspond to models M005V0 (a = 0), MOOS VI (a = IQ-^) and M005V2 (a = 3 x IQ-^), 
respectively. Viscosity makes it easier for gas flows towards the protoplanet. Therefore, the surface 
density of the circumplanetary disk in the viscous gas is higher than that of the inviscid gas. We have 
calculated the circumplanetary disk mass. For a 8.5Mq protoplanet, when a = 3 x 10^'^, the disk mass 
is bigger than the inviscid disk by a factor of 56%. For a 16. 5^/^ protoplanet, when a = 3 x 10"'^, the 
disk mass is bigger than the inviscid disk by an alterative factor of 50%. We find that when Mp > 33M0, 
the effect of viscosity on circumplanetary disk mass is negligibly small. 

3.4 Mass accretion rate 

We talk about the mass accretion rate in this section. We talk about mass accretion rate in real unit (equa- 
tions (11)-(15)). We assume that the planets locate at Qp = 5.2AU. Also, we assume that Mc = ^Mq 
and L = IL©. In this paper, we find that the mass accretion rates for non-viscous models M0026V0 
(8.5M®) and M005V0 (16.5Afe) are 2.2 x 10~^Mjup/yr and 5.0 x lQ~^Mjup/yr, respectively. 
Almost all of the numerical settings of Tanigawa & Watanabe (2002) are same as that in this paper. The 
only difference is that their is no viscosity in Tanigawa & Watanabe (2002). In Tanigawa & Watanabe 
(2002), the smallest mass used is 22Mq. Therefore, we can not directly compare our result to theirs. 
Fortunately, in equation (19) in Tanigawa & Watanabe (2002), they give the dependence of mass accre- 
tion rate on the planet mass. Using their equation (19), the mass accretion rate for 8.5M0 and 16.5Mq 
are 2.02 x 10^^ Mjup/yr and 4.8 x 10^^ Mjup/yr, respectively. Our results for non-viscous gas are 
consistent with that in Tanigawa & Watanabe (2002). 

If viscosity is included, we find that the mass accretion rate becomes higher. The mass accretion rate 
for the viscous models M0026V2 (8.5M© and a = 3 x IQ-^) and MOOS V2 (16.5Afe and a = 3 x IQ-^) 
are 4 x lO^^Mjup/yr and 8.0 x 10~^Mjup/yr, respectively. D'Angelo et al. (2002) using global 
simulations to study the gas flow onto planets. For a 6.4A/0 planet, when a = 4 x 10"'^, they find 
that the mass accretion rate is 1.2 x lO^^Mjup/yr (Fig. 25 in their papers). For a 15M0 planet, when 
a = 4 X lO^'^, they find that the mass accretion rate is 1.5 x 10^^ Mjup/yr (Fig. 25 in their papers). 
Given that the viscosity and planet mass are comparable, it seems that the accretion rate found in this 
paper is much higher than that obtained in D' Angelo et al. (2002). The reason may be due to the 
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depletion of the gas inside the orbital radius of the planet in the global simulations in D' Angelo et al. 
(2002). 

The Bondi accretion rates for 8.5M® and 16. SM® planets located at 5.2AU are 1 x 10~'^Mjup/yr 
and 3.8 x lO^^Mjup/yr, respectively. The angular momentum of the gas makes the actual accretion 
rate much smaller than the Bondi accretion rate of the planets. We must note that the mass accretion rate 
obtained in this paper may be not accurate. This is because, inside the Hill radius, the gas may evolve 
adiabatically, the increased temperature towards the planet may make accretion rates smaller than that 
obtained in this paper 

4 SUMMARY AND DISCUSSION 

We investigate the effects of viscosity on the circumplanetary disks forming in the vicinity of the pro- 
toplanet through two-dimensional hydrodynamical simulations with the shearing sheet model. We find 
that viscosity can affect the properties of the circumplanetary disk significantly when the mass of the 
protoplanet < 33 A/®. 

Local shearing-sheet approximation is only an approximation of the global model, and it may not be 
appropriate for investigating the global evolution of the disk structure. However, Muto et al. (2010) have 
shown that the local shearing-sheet approximation and full global model share many essential physics in 
common. Muto et al. (2010) have also shown that the one-dimensional disk evolution model constructed 
from the global model and the local model are very similar Therefore we expect the local simulations 
have captured the main physics of the circumplanetary disks. 

Physically we should use three-dimensional simulations to study the circumplanetary disks inside 
the Hill radius. However, as a first step, we carry out these simulations in order to understand the basic 
effects of viscosity on the circumplanetary disks. The two-dimensional simulations can capture the 
main physics of the disks. For example, the mass accretion rate can be generated properly by two- 
dimensional simulations because the accretion rate is determined mainly by the flow where r > Rh 
(when r > Rh, the disk is very thin and two-dimensional simulations is enough) (Tanigawa & Watanabe 
2002). Although the 2D simulations can capture some properties of the circumplanetary disk, some 
important feature may lost in the 2D disks. For example, Tanigawa, Ohtsuki & Machida (2012) found 
that most of gas accretion onto circumplanetary disks occurs nearly vertically toward the disk surface 
from high altitude, which can not be found in 2D simulations. In a subsequent paper, we will discuss 
the effects of viscosity on the circumplanetary disks in three-dimensions. 

In this paper, we use a smoothing length to smooth the gravity close to the planet. Muller (et al. 
2012) find that for longer distances, the smooth length is determined solely by the vertical disk thickness. 
For the planet case they find that outside r ^ H, when the value of Vsm = 0-7H describes the averaged 
force very well. However, for shorter distances the smoothing needs to be reduced significantly. In this 
paper, in order to study the structure of circumplanetary disks we adopt rsm ~ Q.05H. This value is 
proved to be safe to study the circumplanetary disks by Tanigawa & Watanabe (2002). 

In this paper, we just use an anomalous stress tensor to mimic the shear stress, which is in reality 
magnetic stress associated with magnetohydrodynamic (MHD) turbulence driven by the magnetoro- 
tational instability (MRI). In a real turbulent circumplanetary disk, the properties of the disk should 
fluctuate with time, but we expect that the time-averaged properties of the turbulent disk should be con- 
sistent with our results here. The spiral shocks in the circumplanetary disk also can affect the amplitude 
of the turbulent stress, but the effect is small when the mass of the protoplanet is small (Mp < SOM^) 
(Papaloizou et al. 2004). Thus, our calculations have captured the main physics of the circumplanetary 
disk. 

To properly study the viscous circumplanetary disk, it is necessary to include magnetic effects. 
However, the gas in the circumstellar disk surrounding a protostar is just weakly ionized. Weakly ion- 
ized plasma is subject to a number of non-idea MHD effects due to the collisional coupling between 
the ionized species and the neutrals (e.g. ambipolar diffusion effects) (Bai & Stone 2011). The mag- 
netorotational instability (MRI), which is considered as the major mechanism for angular momentum 
transport via the MHD turbulence, is strongly affected by the non-idea MHD effects. Thus, it is neces- 
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sary to study the non-idea MHD effects on the MRI before further investigating the properties of the 
magnetized circumplanetary disks. 

We find that when Mp < 33 Af^, viscosity can disrupt the spiral structure around a protoplanet 
considerably and make the gas in the disk smoothly distributed, which makes the torques exerted on 
the protoplanet weaker Thus, viscosity can make the migration speed of a protoplanet lower This is 
helpful to solve the problem that a protoplanet quickly migrates to the vicinity of the central star before 
becoming a gas giant planet. 

According to the core accretion theory, the formation process of a gas giant planet can be divided 
into three phases. (1) In phase 1, the solid core forms first, which has a mass of several Earth mass. 
(2) In phase 2, a spherically hydrostatic gas envelope around the solid core forms. The gas accretion 
rate in this phase is very low. (3) After the point when the core mass and the envelope mass become 
comparable, gas is accreted in a runaway fashion. The main problem of the core accretion theory is 
that the formation time of a gas giant planet exceeds the lifetime of the circumstellar disk. Under the 
assumption of spherical symmetry and gas in hydrostatic equilibrium, previous works found that the 
time needed to complete phase 2 is comparable or exceeding the lifetime of the circumstellar disk (e.g. 
Pollack et al. 1996). The reason for the longtime evolution in phase 2 is that they assume spherical 
accretion. The gravitational energy released in the accretion process can not easily escape from the 
system, the thermal pressure can support the envelope against the gravity of the core, which decreases 
the accretion rate in phase 2 significantly. Lin (Lin 2006) proposed that if a circumplanetary disk (instead 
of a spherically hydrostatic envelope) exists around a several Earth mass protoplanet, the evolution time 
of phase 2 may be decreased significantly. This is because in the disk accretion case, the gravitational 
energy released in the accretion process can easily escape from the surface of the circumplanetary disk. 
However, our simulation found that when the protoplanet is small (several Earth mass), the size of the 
circumplanetary disk is just ~ 27% of the Hill radius, Bh, when viscosity is included (see Fig. 4). Thus, 
even though the existence of a circumplanetary disk can make the formation time of a gas giant planet 
smaller, we expect that the effect is not important because the size of the disk is small enough compared 
to the Hill radius. 

Satellites may form in the circumplanetary disk. Regular satellites around a gas giant planet are 
considered to form according to a scenario similar to the formation of Earth-like planets in our Solar 
system. The protosatellite forms by the accumulation through mutual collision of the satellitesimals that 
form after the dust grains sink towards the equatorial plane of the circumplanetary disk (e.g. Stevenson, 
Harris & Lunine 1986). After including viscosity, the density of the circumplanetary disk increases, 
which should be helpful for the formation of satellitesimals by the accumulation of dust grains. However, 
the increase in density will result in the increase of opacity. It is difficult for high opacity circumplanetary 
disks to radiate its energy generated by gas accretion. Thus, the increase of density may result in the 
increase in temperature of the circumplanetary disks. Large fraction of ice is found in some of the 
Galilean moons around the Jupiter (Schubert, Spohn & Reynolds 1986). This means that the temperature 
of the circumplanetary disk should be low enough for water condensation into ice (Canup & Ward 2002). 
In this sense, the increase of density may make it difficult for satellite formation. In order to understand 
the properties of the circumplanetary disk better, we need to study the viscous circumplanetary disk with 
radiative transfer, which is beyond the scope of this paper. 
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